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The voids between galaxies are identified with the volumes of the Poisson 
Voronoi tessellation. Two new survival functions for the apparent radii of 
voids are derived. The sectional normalized area of the Poisson Voronoi tes- 
sellation is modelled by the Kiang function and by the exponential function. 
Two new survival functions with equivalent sectional radius are therefore de- 
rived; they represent an alternative to the survival function of voids between 
galaxies as given by the self-similar distribution. The spatial appearance of 
slices of the 2dF Galaxy Redshift Survey is simulated, 
keywords 

methods: statistical ; cosmology: observations; (cosmology): large-scale 
structure of the Universe 



1 Introduction 

The statistical analysis of the voids between galaxies is a field of research 
and the following catalogs have been explored: the two-degree Field Galaxy 
Redshift Survey (2dFGRS), see Patiri et al. (2006) and von Benda-Beckmann 
& Miiller (2008), the Sloan Digital Sky Survey (SDSS), see Tikhonov (2007), 
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and the combination of 2dFGRS and SDSS, see Tinker et al. (2008). The 
voids between galaxies are usually modelled by a survival function in the ap- 
parent radius as given by a modified exponential distribution, see (3) in von 
Benda-Bcckmann & Miiller (2008) or our Section 4; this fact is considered 
here a standard argument for comparison. The concept of Voronoi Diagrams 
dates back to the vortex theory applied to the solar system as developed in 
the 17th century, see Descartes (1644) and Figure 1 in Aurenhammer, F. 
and Klein, R. In : Sack, J.-R. (2000). The name is due to the two historical 
records by Voronoi (1907, 1908) and the applications to astrophysics begin- 
ning with Kiang (1966). The Voronoi diagrams represent a model of the 
voids between galaxies, sec van de Wcygaert & Icke (1989), Pierre (1990), 
Barrow & Coles (1990), Coles (1991), van de Weygaert (1991), van de Wey- 
gaert (1991), Zaninetti (1991), Ikeuchi & Turner (1991), Subba Rao & Szalay 
(1992), van de Weygaert (1994), Goldwirth et al. (1995), van de Weygaert 
(2002), van de Weygaert (2003), and Zaninetti (2006). 

The Poisson Voronoi tessellation (PVT) is a particular case of the Voronoi 
tessellation in which the seeds are generated independently on the X, Y and 
Z axes in 3D through a subroutine which returns a random real number 
taken from a uniform distribution between and 1. For practical purposes, 
the subroutine RAN2 was used, see Press et al. (1992). On adopting an 
astrophysical point of view, the sectional PVT, Vp{2,3), is very interesting 
because it allows a comparison with the voids as observed in slices of galaxies 
belonging to different catalogs such as the CFA2 catalog (Geller & Huchra 
(1989)), the 2dfGRS (Norberg et al. (2002)), the 6dF Galaxy Survey (GdFGS) 
(Jones et al. (2004)) or the SDSS (Einasto et al. (2003)). The absence of (i) 
a numerical analysis through the survival function of normalized areas in 
2D and normalized volumes in 3D of PVT (ii) a careful exploration of the 
statistical properties of \^(2, 3), leads to the following questions being posed. 

• Is it possible to integrate the usual probability density functions (PDFs) 
which characterize the main parameters of 2D and 3D PVT in order to 
obtain an analytical expression for the survival function? 

• Is it possible to model the normalized areas of Vp{2, 3) with the known 
PDFs? 

• Can we transform the normalized volumes and areas into equivalent 
radius distributions? 



2 



• Can we simulate the observed slices of galaxies as given, for example, 
by the 2dF Galaxy Redshift Survey? 

In order to answer these questions. Section 2 reports the three major 
PDFs used to model the normahzed area/volume of 2D/3D PVT as well as 
the results of the fit. 

Section 3 reports the apparent distribution in effective radius of the 3D 
PVT as well as their associated survival functions. 

Section 4 contains the self-similar survival function for voids between 
galaxies as well as the associated PDF. 

Section 5 reports the fit of the normalized area distribution of the sectional 
PVT with the Kiang function and the exponential distribution. 

Section 6 reports our actual knowledge of the photometric properties of 
galaxies as well as a Voronoi simulation. 

It is important to point out that the PVT is not used as a technique to 
identify voids in existing data catalogs, see Ebehng & Wiedenmann (1993), 
Bernardeau & van dc Wcygaert (1996), Schaap & van dc Wcygaert (2000), 
Marinoni et al. (2002), Melnyk et al. (2006), van de Weygaert & Schaap 
(2009) and Elyiv et al. (2009). 

The PVT formalism is here used conversely: to generate mock catalogs 
which are later calibrated by observations. On adopting the point of view 
of the statistical distributions is important to underline that the survival 
fimction is here identified with the cumulative void size distribution function. 
We briefiy recall that the cumulative void size distribution function relates 
the number of voids to their size, analoguosly to the halo mass function which 
relates number of dark matter halos to their mass. 



2 The adopted distributions of the PVT 



A PDF is the first derivative of a distribution function (DF) F{x) with respect 
to X. In the case where the PDF is known but the DF is unknown, the 
following integral is evaluated 




(1) 



As a consequence the survival function (SF) is 

SF^l- F{x) 



(2) 
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2.1 The Kiang function 

The gamma variate H{x; c) (Kiang (1966)) is 

H{x;c) = Y^{cxy'^exp{-cx) , (3) 

where < x < oo, c > 0, and r(c) is the gamma function with argument c. 
The Kiang PDF has a mean of 

/^=1 , (4) 

and variance 

.-1 . (5) 

In the case of a ID PVT, c = 2 is an exact analytical result and conversely 
c is supposed to be 4 or 6 for 2D or 3D PVTs, respectively, Kiang (1966). 
The DF of the Kiang function, DFk, is 

^ . (a) 

where the incomplete Gamma function, r{a,z), is defined by 

r(a, z) = / e-Y~^dt . (7) 
The survival function Sk is 

r (c, cx) 

2.2 Generalized gamma 

The generalized gamma PDF with three parameters a, b, c, Hinde & Miles 
(1980); Ferenc & Neda (2007); Tanemura (2003), is 

f^a/c 

f{x;b,c,d)^c x" 'expi-bx") . (9) 

r(a/c) 



4 



The generalized gamma has a mean of 



/^=^77^ > (10) 



and a variance of 

, ^ ^-H+r(^)r(f)-(r(?;, ^ 

The DF of the generahzed gamma is 

DFgg^I- J . (12) 

The SF of the generahzed gamma is 

ft 



Sgg= 1%^^ . (13) 



2.3 Ferenc— Neda function 

A new PDF has been recently introduced, Ferenc & Neda (2007), in order 
to model the normalized area/volume in a 2D/3D PVT 

FN{x;d) = C X x'^ exp{-{3d+l)x/2) , (14) 
where C is a constant, 

C = ^l^gH . (15) 

2 23/2 (3 d + 1 ) r (3/2 + 1/2) 

and d{d = 1, 2, 3) is the dimension of the space under consideration. We will 
call this function the Ferenc-Neda PDF; it has a mean of 

/x = l , (16) 
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and variance 

2 



3d+l 

The DF of the Ferenc-Neda function when d — 2 is 

DF, 



FN2 



49 V2V7x^/^e-y^^ ^/2V7x^/^e-y^^ 
15 

-^?^I|^ + erf(l/2V2V7V^), 
where the error function erf (a;) is defined as 

erf (a;) = / 2 -^dt . 

Jo VTT 



The SF of the Ferenc-Neda function when d = 2 is 



S 



FN2 — 

1 + — ^-^^ T= + 7/3 ^ 

15 VTT VTT 



The DF of the Ferenc-Neda function when o? = 3 is 



DFpm — 

1 - e"^^ - 5e~^''a; - —e'^'-'x^ 

2 

125 3 625 4 

e X X e . 



6 24 
The SF of the Ferenc-Neda function when c? = 3 is 



+^e-^V + ^e-5-x^ + ^a;V^- . 
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Table 1: Values of for the cell normalized area-distribution function in 
2D; here Tj are the theoretical frequencies and Oj are the sample frequencies. 
Here, we have 25087 Poissonian seeds and 40 intervals in the histogram. 



PDF 


parameters 


x' 


H{x;c){Eq.{S)) 


c = 3.55 


83.48 


f{x;d){Eq.{U)) 


d = 2 


71.83 


G{x;a,b,c){Eq.{9)) 


a = 3.15 


58.9 




6 = 2.72 






c = 1.13 





2.4 Numerical results 

In the following, we will model the PVT in which the seeds are computed 
through a random process. The is computed according to the formula 

i=l 

where N is the number of bins, Tj the theoretical value and Oi the experi- 
mental value. A first test of the PDFs presented in the previous section can 
be done by analysing the Voronoi cell normalized area-distribution in 2D, see 
Table 1. 

Table 2 reports the parameters of the PDF of the volume distribution in 

3D. 

Figure 1 reports an example of SF in 2D PVT (areas) and Figure 2 an 
example of SF in 3D PVT (volumes) . 

3 Radius distribution of the 3D PVT 

We now analyse the distribution in effective radius R of the 3D PVT. We 
assume that the volume of each cell, v, is 

v-U-r , (24) 
3 p 
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Table 2: Values of foi^ cell normalized volume-distribution function in 3D. 
Here, are the theoretical frequencies and Oi are the sample frequencies. 
Here, we have 21378 Poissonian seeds and 40 intervals in the histogram. 



PDF 


parameters 


x' 


H{x-c) {Eq.{3)) 


c = 5.53 


93.86 


f{x-d) {Eq.{U)) 


d = 3 


134.15 


G{x;a,b,c) iEq.{9)) 


a = 4.68 


58.59 




6 = 3.87 






c= 1.18 






Figure 1: SF of normalized area-distribution function in 2D when we have 
25087 Poissonian seeds and 40 intervals: the empty circles represent the 
Voronoi volumes and the full hne the theoretical value of Sgg (generalized 
gamma function) as represented by (13). 
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Figure 2: SF of normalized volume-distribution function in 3D when we 
have 21378 Poissonian seeds and 40 intervals: the empty circles represent 
the Voronoi volumes and the full line the theoretical value of Sk (Kiang 
function) as represented by (8). 



where p is a length that connects the normalized radius to the observed one. 
In the following, we derive the PDF for radius and related quantities relative 
to the Kiang function and Ferenc-Neda function. 

3.1 Kiang function of the radius 

The PDF, Hr{R;c), of the radius corresponding to the Kiang function as 
represented by (3) is 



Hn{R;c) = ^ -YW? ' ^^^^ 

where < R < oo, c >0 and p > 0. The Kiang PDF of the radius has a 
mean of 

^ ^^r(i/3 + c) 
^" / ^^r(c) ^ ' ^^^^ 

and variance 
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1 3i2§ (r (2/3 + c) r (c) - (r (1/3 + c)f) 



4 c2/37r2/3(r(c)) 
The DF of the Kiang function of the radius is 

_^ r(c,4/3cvr(y) 



r c 



The survival function of the radius is 

r(c,4/3c7r(ff 



Skr — 



r(c) 

3.2 The Ferenc— Neda function of the radius 

The PDF as a function of the radius, obtained from (14) and inserting 
is 

400000 TT^i^^^e" 



^ 1 . _20 

-5pl4- 3 ^ 



FNr{R; d) 



243 

The mean of the Ferenc-Neda function is 

IX = 0.6p , 

and the variance is 

= 0.0085p2 . 

The DF of the Ferenc-Neda function of the radius when d = 3 is 

DFpNSR = 

_2() TTfl^ 20 -20 7rfl3 

1 - e ' - y e 3 ~^R^TTp-^ 
200 . 2 fi 
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81 
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243 ^ 
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The SF of the Ferenc-Neda function of the radius when c? = 3 is 

SfN3R — 

_ 20 7rfli 20 -20 TTfli 

e 3 _^ g 3 
9 

4000 _2o^ 

ol 

20000 -20^ 

H e ^ ^i?^Vp-^2 . (34 

243 ^ ^ ^ 

4 Self-similar void distribution 

The statistics of the voids between galaxies have been analysed in von Benda- 
Beckmann & Miiller (2008) with the following self-similar SF in the following, 

Sss, 

Sss^e V^iV , (35) 

where A is the mean separation between galaxies, si and S2 are two length 
factors, and pi and p2 two powers. The DF of the self-similar distribution is 



DFss = 1 - c V^iV V=2V _ (36) 
The PDF of the self-similar distribution is 

Pss{R) = 

R 

At present, it is not possible to find an analytical expression for the integral 
that defines the average value of the self-similar distribution. 

A comparison of the survival function of self-similar voids with the sur- 
vival function of the radii of the two PDFs analysed here is reported in 
Figure 3. 

The two fits of Figure 3 are not satisfactory because we are making a 
comparison between the radius of projected voids and the 3D radii of the 
normalized volume distribution of the PVT. This fact is confirmed from the 
high values of computed according to (23). 
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Figure 3: The survival function, Sss, of the self-similar distribution in radius 
of N/Sl as represented by (35) (full hue), the survival function, Sfnsr, for 
the Ferenc-Neda function of the radius as represented by (34) when d — 3, 
p — 12.5 Mpc and = 4319 with 100 points (dashed line). The survival 
function, Skr, of the radius of the Kiang function as represented by (29) 
when p = 12.5 Mpc, c = 5.53 and = 4076 with 100 points (dot-dash-dot- 
dash line). 
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Tabic 3: Values of for the cell normalized area-distribution function of 
Vp{2,3); here Tj are the theoretical frequencies and Oj are the sample fre- 
quencies. Here we have 8517 Poissonian seeds and 40 intervals in the his- 
togram. 



PDF 


parameters 




H{x;c) {Eq.{3)) 


c = 2.07 


114.41 


p{x;b) {Eq.{38)) 




85.38 



Table 4: Average value of c according to the Kiang function, Eq. (3)), for 
Vp{2, 3), when 50 cuts are considered. Here we have 25087 Poissonian seeds. 

PDF name c a 

H{x;c) {Eq.{3)) Kiang 2.25 0.11 



5 The sectional PVT 

The existence of voids between galaxies is normally deduced from a projected 
distribution in astronomical slices such as the 2dFGRS S3, see Figure 1 in 
von Benda-Beckmann & Miiller (2008). This fact motivates the analysis of 
the sectional Voronoi tessellation, also known as Vp{2,3), see Okabe et al. 
(1992). A previous analysis has shown that a cell belonging to the intersec- 
tion between an arbitrary plane and the faces of the Voronoi polyhedrons 
is almost surely a non- Voronoi cell, see details in Chiu et al. (1996). Here, 
we first model the normahzed area-distribution 1^(2, 3) with Kiang PDFs as 
represented by (3), see Table 3. 

In the case of n cuts, we can compute the average value of c, c, and the 
standard deviation a, see Table 4. 

A test of this value can be done on the unpublished manuscript of Ken 
Brakke available at http://www.susqu.edu/brakke/aux/ 
downloads/papers/3d.pdf . 

Table 4 of this paper is dedicated to the 3D plane cross-sectional statistics 
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E(area) = 0.6858 and Var(arca) = 0.2269. When a typical run of ours is 
normalized to the average value rather than 1, our result is Var(area) = 
0.2266 which means that our numerical evaluation differs by 2/10000 from 
the theoretical result. We remember that in the case of the normalized area- 
distribution function in 2D we have c ~ 4 and we can therefore speak of a 
decrease in the value of c by a factor 2. The decrease of c for Vp{2, 3) was 
first derived in Figure 4 of Zaninetti (2006). 

Another PDF that can be considered in order to model the normalized 
area distribution of Vp{2, 3) is the exponential distribution, 

p{x) = ^exp~ , (38) 

which has an average value 

x^b . (39) 

In the case of the normalized areas b = 1, Table 3 reports the values of the 
two distributions adopted here. Once the statistics of the normalized area 
distribution of an arbitrary cut Vp{2, 3) are known, we can model the radius 
distribution. We therefore convert the area of each cell. A, to an equivalent 
radius R 

R=\ - . (40) 
V vr 

We also briefly remember that the problem of stereology is to deduce the true 
size distribution N{R) of 3D-volumes from the distribution n(r) of apparent 
circles in 2D. 



5.1 Kiang distribution of 1^(2,3) in radius 

The PDF, Hr2z{R] c), as a function of the radius corresponding to the Kiang 
function as represented by (3) for V^(2, 3) is 

2c(^)'^-^-"^7ri? 
Hm^{R;c)= ^ ' I.. , , (41) 

r (c) 

where 0<i?<oo,c >0 and p > 0. The Kiang PDF of the radius for 
Vp{2, 3) has a mean of 

^-n^ . (42) 

^fc^jT^Y (c) 
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and variance 

, ^ p'{c{T{c)f-{T{c+l/2)f) 

C7r(r(c))^ ■ ^ ^ 

The DF of the Kiang function of the radius, DFkr23, for Vp{2, 3) is: 

The survival function Skr22 of the radius for V^(2, 3) is 

r (c, 2 ^) 



r(c) 



(45) 



A comparison of the survival function of self-similar voids with the sur- 
vival function of the radius for V^(2, 3) of the exponential distribution is 
reported in Figure 4. 

5.2 Exponential distribution of Vp{2, 3) in radius 

The PDF, Pr2z{R] c), as a function of the radius corresponding to the expo- 
nential distribution as represented by (38) for Vp{2, 3) is 

,^ , 2e ~p^ttR 

PR23{R;c)^ , (46) 

p 

where < i? < oo, p > 0. The exponential PDF of the radius for 1^(2, 3) 
has a mean of 

7^3/2 

. .3/2 ' (47) 

and variance 

(48) 



An 

The DF of the exponential distribution in radius, DFer23: for 1^(2, 3) is: 

,2 



DFER23^1-e ^ . (49) 
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Figure 4: The survival function, Sss, of the self-similar distribution in radius 
of N/Sl as represented by (35) (full line), the survival function, Skr23: of 
the radius of the Kiang function for 1^(2, 3) as represented by (45) when 
p = 13 Mpc, c = 2.25 and = 67.1 with 100 points (dot-dash-dot-dash line). 
The survival function, Ser23, of the radius of the exponential distribution for 
Vp{2, 3) as represented by (50) when p = 7 Mpc and = 9-27 with 100 points 
(dashed line). 
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The survival function Ser23 of the radius for Vp{2, 3) is 

Ser23 = e . (50) 

Figure 4 reports a comparison between the survival function of self-similar 
voids and the exponential distribution for 1^(2,3). 

In this case, the two types of fit of Figure 4 are satisfactory because 
we are making a comparison between the observed projected radius of voids 
and the projected radius of the Voronoi volumes. The smaller associated 
with the exponential distribution can be considered a consequence of the fact 
that the statistics of the normalized area distribution of the cuts are better 
described by an exponential distribution than by the Kiang function, see 
Table 3. 

A final comparison between the four samples of void size statistics as 
represented in Figure 4 of von Benda-Beckmann & Miiller (2008) and our 
survival function of the radius of the exponential distribution for Vp{2, 3) is 
reported in Figure 5. 

6 Large-scale structure 

We briefly review the status of the knowledge of the Hubble constant, refer- 
ence magnitude of the sun, luminosity function of galaxies, Malmquist bias 
and 3D Poissonian Voronoi diagrams. 

6.1 The adopted constants 

A first important evaluation of the Hubble constant is through Ccphcids 
(key programs with HST) and type la Supernovae, see Sandage et al. (2006), 
Hq = (62.3 ± 5) km s~^ Mpc~^. A second important evaluation comes from 
three years of observations with the Wilkinson Microwave Anisotropy Probe, 
see Table 2 in Spergel et al. (2007), Hq = (73.2 ± 3.2) km s~^ Mpc~\ In the 
following, we will take the average value of these two important evaluations, 
Hq = 67.65 km s~^ Mpc~^. The Hubble constant is also reported as 
Hq = lOO/i km s~^ Mpc~^, with h = 1 when h is not specified; in our case 
h = .6765. 
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Figure 5: The survival function, Sss, for the self-similar distribution in radius 
of N/S2, N/S3 and as reported in Figure 4 of von Benda- 

Beckmann and MuUer 2008 (full red lines), as represented by (35). The 
survival function, Ser23, of the radius of the exponential distribution for 
Vp{2, 3) as represented by (50) when p = 7 Mpc, p = 9 Mpc, p = 13 Mpc 
and p = 17.7 Mpc (dashed green lines). 
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Another quantity that should be fixed in order to continue is the absolute 
magnitude of the sun in the bj filter of the 2dFGRS M.Q = 5.33, see CoUess 
et al. (2001); Tempel et al. (2009); Eke et al. (2004). 

6.2 Malmquist bias 

This bias was originally applied to the stars, see Malmquist (1920, 1922), 
and was then applied to the galaxies by Behr (1951). We now introduce the 
concept of limiting apparent magnitude and the corresponding completeness 
in absolute magnitude of the considered catalog as a function of redshift. 
The observable absolute magnitude as a function of the limiting apparent 
magnitude, m^, is 

Mi = mL-51og,o(^)-25 . (51) 

The previous formula predicts, from a theoretical point of view, an upper 
limit on the absolute maximum magnitude that can be observed in a catalog 
of galaxies characterized by a given limiting magnitude. 

The interval covered by the LF of galaxies, AM, is defined by 

AM = M^„^-M^i„ , (52) 

where M^nax ^-nd M^j„ are the maximum and minimum absolute magnitude 
of the LF for the considered catalog. The real observable interval in absolute 
magnitude, AMl, is 

AML = ML-M^i^ . (53) 

We can therefore introduce the range of observable absolute maximum mag- 
nitude expressed in percent, e{z), as 

^.(^) = ^xlOO% . (54) 

This is a number that represents the completeness of the sample and, given 
the fact that the limiting magnitude of the 2dFGRS is 171^=19.61, it is possi- 
ble to conclude that the 2dFGRS is complete for z < 0.0442 . This efficiency, 
expressed as a percentage, can be considered to be a version of the Malmquist 
bias. 
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Table 5: Parameters of the Schechter function forthe 2dFGRS. 



parameter 2dFGRS 

M* - 5 logio h [mags] (-19.79 ± 0.04) 

a -1.19±0.01 

$* [h^ Mpc-^] ((1.59 ± 0.1)10-2) 



6.3 Luminosity function of galaxies 

The Schechter function, introduced by Schechter (1976), provides a useful fit 
for the LF (luminosity function) of galaxies 

mdL = {^){^reM-^)dL . (55) 

Here, a sets the slope for low values of L, L* is the characteristic lumi- 
nosity and $* is the normalization. The equivalent distribution in absolute 
magnitude is 

^{M)dM = (0.4ZnlO)$nO°-^("+i)(^*-^) 

xexp(-10°-^(^*-^))dM , (56) 

where M* is the characteristic magnitude as derived from the data. The 
parameters of the Schechter function for the 2dFGRS can be found on the 
first line of Table 3 in Madgwick et al. (2002) and are reported in Table 5. 

We now introduce /, the flux of radiation expressed in units of Lq Mpc"^, 
with Lq representing the luminosity of the sun. The joint distribution in z, 
redshift and /, see (1.104) in Padmanabhan (1996) or (1.117) in Padmanab- 
han (2002), is 

where dfl, dz, and df represent the differential of the solid angle, redshift, 
flux and q represents the velocity of light. 

This relationship has been derived assuming z ^ ^ with r rep- 
resenting the distance of the galaxy in Mpc. The critical value of z, Zcrit, 
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IS 



^'"*"4^7^ ■ ^^^^ 

The number of galaxies in z and / as given by (57) has a maximum at 

z — ^pos—max ^ where 

= Zcrit^/Oi^ 2 , (59) 



z, 



pos—max 



which can be re-expressed as 

'^pos—max — r> / — nn ■ V / 

2 V^VJ Q 

The number of galaxies, Ns{z, fmin, fmax), comprised between the min- 
imum value of flux, fmin, and the maximum value of flux, fmax, can be 
computed with the following integral 



2 



fmax fir. 7 

^s{z) = 47r(^)V$(^)rf/ . (61) 



fmin Hq Z^ 



2 

crit 



This integral does not have an analytical solution and therefore a numerical 
integration must be performed. The total number of galaxies in the 2dFGRS 
is reported in Figure 6 as well as the theoretical curves as represented by 
the numerical integration of (57). A careful inspection of the previous figure 
allows to conclude that the theoretical integral fits well the experimental 
data up to 2; = 0.0442. Beyond that value, the presence of the Malmquist 
bias decreases the number of observable galaxies and the comparison between 
theory and observations cannot be made. 



6.4 Voronoi diagrams 

Voronoi diagrams represent a useful tool for describing the spatial distribu- 
tion of galaxies and, as an example, van de Weygaert & Icke (1989) identified 
the vertexes of irregular Voronoi polyhedrons with Abell clusters. Another 
example is provided by Zaninetti (1991) where the galaxies were first inserted 
on the faces of the irregular Voronoi polyhedrons and a power law distribu- 
tion for the seeds was adopted. Later, the galaxies were still inserted on the 
faces of the irregular Voronoi polyhedrons but Poissonian seeds were adopted, 
see Zaninetti (2006, 2008, 2010). Adopting the same algorithm developed in 
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Figure 6: The galaxies in the 2dFGRS with 13.34 < bJmag < 19.61 or 
1540 Lq Mpc~^ < / < 493844L0 Mpc"^, are organized as frequencies 
versus hehocentric redshift (empty stars). The theoretical curve generated 
by the integral of the Schechter function in flux, (61, (full line) is drawn. 
The maximum in the frequencies of observed galaxies is at 2; = 0.108. In 
this plot, A^0 = 5.33 and h = 0.67. The vertical dotted line represents the 
boundary between complete and incomplete samples. 
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Mpc 



1000 



Figure 7: A Voronoi slice 75° long and 3° wide. The range of the flux is 2500 
Lq Mpc-2 < / < 2868O8L0 Mpc'^ and the number of galaxies is 19895. 
The red circle denotes the confusion distance as given by (64). 
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Zaninetti (2010), one slice of the 2dFGRS with the number of galaxies as a 
function of z as given by (61) is simulated, see Figure 7. 

This simulation should be compared with Figure 1 in von Benda-Beckmann 
& Miiller (2008) and with a 3° slice of the 2dFGRS Image Gallery available 
at the web site, http://msowww.anu.edu.au/2dFGRS/. Little need be said 
about the number of seeds which should be used in order to simulate the 
observed slices of the 2dFGRS. The average volume of the voids is the side 
of the box in Mpc divided by the number of seeds. The average diameter of 
the voids, Dy^anm-, is 

^^'-""''^ - of seeds^ " ^^^^ 

The theoretical average diameter, D^h, can be obtained from the average 
value (47 ) of the exponential distribution of Vp{2, 3) in radius once the 
maximum value of p which fits the von Benda-Beckmann & Miiller (2008) 
data, see Figure 5, is adopted 

Dth^f^l . (63) 

The equality Dvoronoi = Dth allows to obtain the number of Poissonian seeds. 

Particular attention should be paid to the fact that the astronomical shces 
arc not a plane which intersects a Voronoi network. In order to quantify this 
effect, we introduce a confusion distance, D^, as the distance after which the 
altitude of the slices equalizes the observed average diameter Dof,s 



DctSin(a) = Do,,s , (64) 



where a is the opening angle of the slice and Dobs is the average diameter 
of the voids. In the case of 2dFGRS a = 3° and therefore = 498.5 Mpc 
when Dobs — 26.12 Mpc, see the red circle in Figure 7. 



7 Conclusions 

The PDFs which are usually used to model the normalized volume distribu- 
tion of the 3D PVT are gamma type distributions such as the Kiang function, 
(3), and the new Ferenc-Neda function, (14). 
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Here, in order to make a comparison with the self-similar distribution of 
voids, we derived the survival distribution of the Kiang function, (8), and of 
the Ferenc-Neda function, (22). 

On adopting an astrophysical point of view the cut Vp(2, 3) may model 
the voids between galaxies as given by astronomical slices of the Millenium 
catalog. The analysis of the normalized areas of V.p{2. 3) is a subject of re- 
search rather than a well-established fact and we have fitted them with the 
Kiang function and the exponential distribution. The value indicates that 
the exponential distribution fits more closely the normalized area distribu- 
tion of Vp{2,3) than does the Kiang function, see Table 3. This fact follows 
from the comparison between the self-similar survival function and the ex- 
ponential and Kiang distributions of the radius, see Figure 4. Therefore, the 
one parameter survival function of the radius of the exponential distribution 
for Vp{2,3), Ser23, as represented by (50), may model the voids between 
galaxies as well as the five parameter self-similar survival function. 

The observed 2dFGRS slices can be simulated but the behaviour of the 
luminosity function for galaxies and the consequent number of galaxies as a 
function of the redshift should be carefully analysed. 

We are planning, in future projects, 

• To insert the thickness of the faces of PVT and to model the connected 
modification of the survival function. 

• To apply the technique here developed to the real data sets. Millennium 
simulation, for example. The 3D nature of the method and detailed 
density mapping of the voids should be superior to other void identifi- 
cation algorithms as suggested by van de Weygaert & Schaap (2009). 
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